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SUMMARY 


A numerical solution of an elastohydrodynamically lubricated (EHL) line 
contact between two long, rough-surface cylinders that considers the fric- 
tional heating of asperities was obtained. Pressure distribution, tempera- 
ture distribution, film thickness and EHL load for given speeds, lubricant 
properties, material properties of surfaces, and surface roughness parame- 
ters were theoretically solved by simultaneous solution of the elasticity 
equation and the Reynolds equation for two partially lubricated rough sur- 
faces. The pressure due to asperity contact was calculated by assuming a 
Gaussian distribution of surface irregularities. The elastic deformation 
used for film thickness computation was found from the two kinds of pressure 
by plane strain analysis. The temperature rise in the contact zone was cal- 
culated by using the Blok-Jaeger flash temperature model. The effect of 
surface roughness on EHL load for various slide-roll ratios, surface rough- 
ness parameters, surface patterns, and temperature parameters was studied. 

It was found (1) that the maximum temperature rise in most cases occurred in 
the inlet zone and (2) that minimum film thickness decreased and maximum 
temperature increased as the surface roughness was increased. 


INTRODUCTION 

In an earlier report (ref. 1) the effect of the surface roughness of an 
isothermal, elastohydrodynamically lubricated (EHL) line contact was presen- 
ted. In heavily loaded, nonconformal machine elements under elastohydrody- 
namic lubrication there is a possibility of mutual interference between the 
two surfaces. Although some energy is introduced in an elastohydrodynami- 
cally lubricated contact by viscous heating and compression of the lubri- 
cant, the predominant thermal energy source is frictional dissipation due to 
asperity interaction. Frictional heating in the thin-film zone increases 
the temperature and consequently decreases the viscosity of the lubricant. 
This gives rise to a reduction in the film thickness. If the temperature is 
high enough, the lubricating film can fail. Hence the study of the thermal 
effect in concentrated contacts is important. 

Probably Blok (ref. 2) was the first person to use the concept of the 
flash temperature of a solid as it passes through a heat source. Subse- 
quently Jaeger (ref. 3) presented a theory for estimating flash temperature, 
and Archard (ref. 4) gave a procedure for determining this temperature. A 
theoretical study on the thermal stability of such a problem has been made 
by Christensen (ref. 5). He showed that there was a possibility of thermal 
instability because of the reduction in film thickness and the increase in 
the degree of interference of the surfaces. Dyson (ref. 6) discussed vari- 
ous models of the thermal instability of rough-surface contacts. More re- 
cently Patir and Cheng (ref. 7) gave a solution for the inlet region of an 
EHL line-contact problem that considered surface roughness and temperature 
effects. 

Experimental studies on the temperature effect in an EHL contact have 
been performed by Orcutt (ref. 8), Hamilton and Moore (ref. 9), and Kannel 
and Bell (ref. 10). They (refs. 8 to 10) used a platinum wire as the tem- 
perature transducer and measured the temperature of smooth surfaces in line 
contact. Hamilton and Moore (ref. 9) found that the maximum temperature 
rise was always at the inlet zone. Kannel and Bell (ref. 10) used a more 
sophisticated transducer and measured both pressure and temperature. Winer 
and his coworkers (refs. 11 to 13) developed an infrared emission technique 


to measure the temperatures in sliding and rolling contacts for both smooth 
and rough surfaces. The experiments showed that the temperature in the in- 
let zone could be as much as 40° C above ambient temperature. 

The present study estimates the minimum film thickness and the contri- 
bution of energy dissipation and load for an EHL line contact while taking 
into consideration the surface roughness effect. Although the method of 
solution is similar to that of reference 1 it considers the variation of 
viscosity due to temperature rise caused by frictional heating. Viscous 
heating and heating due to the compressibility of the lubricant are neglec- 
ted. The effects of sliding as well as rolling speeds, roughness parame- 
ters, load, lubricant properties, material properties, and the roughness 
patterns of EHL line contacts are investigated. 


SYMBOLS 

B,b 

Ea.Eb 

E' 

G 

H,h 

Hmin» l^min 

Ho»ho 

hf 

K 

k 

P, P 
Pc’Pc 

Q, q 
R 

Pa’^b 

Rs 

T 

To 

Tb 

Tc 

u,ua,ub, 

U , Ug , Ug 

V, v 

Vra’^^rb 
l^c ’ '^c 

Wh,Wh 

Wj,wj 

s,x,y,^ 

S,X,X 

a 

6 

eo» 6 i 

Y 


half the Hertzian contact length, B = b/R, m 
specific heat of metal surface, 0/kg K 
modulus of elasticity of material of surfaces a and b, 
N/m^ 

composite modulus of elasticity, 

1/E' = 1/2 [(1 - v|)/Ea + (1 - vg)/EbJ, N/m2 
materials parameter, G = aE' 
nominal film thickness, H = h/R, m 
minimum film thickness, Hfj,-j,^ = h^-j^/R, m 
central film thickness, Hq = hp/R, m 
average film thickness, m 
a constant 

thermal conductivity of metal, J/m K s 

mean hydrodynamic pressure, P = p/E‘, N/m^ 

contact pressure. Pc = Pc^E', N/m^ 

modified pressure, Q = q/E', N/m'^ 

equivalent radius of cylinder pair, 1 /R = (1/Ra) 

( 1 /Rb)» m 

radii of cylinders a and b, m 
slide-roll ratio, Rs = u^/u 
temperature, K 

temperature of oil at inlet condition, K 

a constant, T 5 = (eo)/(To + ei)^ 

heating parameter, = ^yjUciVR E‘/[BoV'"pck 

X (Vn^ - V%)] ' 

velocity of surfaces, u = (ug + u^)JZ; Uc = 

(ug - U 5 )/ 2 ; U = ngu/E'R; U 5 = t\qu^IE'R, m/s 
elastic deformation of surfaces, = v/R, m 
variance ratios, ( 03 / 0 )^, (ob/o )2 

asperity contact load per unit length, = Wp/E'R, N/m 
hydrodynamic load per unit length, Wp, = w^/E'R, N/m 
total EHL load per unit length, Wj = wj/E'R; 

Wj = W(i Wc, H/m _ 

coordinates, S = s/R; X = x/R; X = x/b, m 

pressure-viscosity coefficient, m^/N 
mean radius of curvature of asperities, m 
constants, K 

surface pattern parameter 
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Ya»Yb 

n^no^ni 

y 

0 

A 

P 

CTa» afa 


‘Py 

<Ps 


THEORY 

For the EHL contact shown in figure 1 the mean hydrodynamic pressure 
and asperity contact pressure are to be determined for various input 
conditions such as loads, speeds, materials, lubricants, and inlet lubricant 
temperatures. 



(a) Two heavily loaded, long rough-surface cylinders. (b) Equivalent cylinder near a plane. 

Figure 1. - Elastohydrodynaniic contact. 


surface pattern parameter of surfaces a and b 
coefficient of absolute viscosity of lubricant, 
no = the same at inlet condition; ni = a constant, N s/m^ 
coefficient of friction between interacting asperities 
dimensionless temperature rise, e = T/pq 
hydrodynamic roughness parameter, a = ho/a 
mass density of metal, kg/m^ 

standard deviations of roughness amplitude of surfaces a and 

b, a = }/ + a§ 

shear flow factor associated with a single surface 
pressure flow factors 

shear flow factor associated with two surfaces 
Poisson's ratio of material of surfaces a and b 
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The mean hydrodynamic pressure is found from the simultaneous solution 
of the average Reynolds equation and the elasticity equation. The governing 
equation for hydrodynamic pressure for rough surfaces when side leakage and 
local squeeze terms are neglected can be written as (ref. 7) 


u + 
a 


u. 3h 


T 


^a - 


ax y**x 

TZri axj " 2 ax 

where 


CPX 

pressure flow factor 

hT 

average gap 

a 

standard deviation 

<Ps 

shear flow factor 

and Hj 

is given by 


3<P; 


( 1 ) 


■X 


(h + 6)f(6)d6 


( 2 ) 


where f(6) is the probability density function of combined roughness 6. 
The nominal film thickness h for this type of contact is approximately 
given by 


|t| _ h + + W 

n - nQ 2R '' 


(3) 


where 




central film thickness 

equivalent radius of curvature of the pair of cylinders 
elastic deformation of cylinders 


The coefficient of dynamic viscosity n is assumed to be varied with 
pressure and temperature according to the following expressions: 


n = Hi e 


ap 


(4) 


and 


m = Tin exp 


J + e- 




( 6 ) 


where 


a pressure-viscosity coefficient 

nQ viscosity at inlet condition where the pressure is atmospheric 

constants 



T To + AT 

Tq inlet temperature 

aT rise in temperature 

By using equation (4) and letting (u^ + ut>)/2 = u and (ug - U[j)/2 = U 5 , 
equation ( 1 ) can be written as 



The pressure flow and shear flow factors 9 ^ 'f’s 

found by Patir and Cheng (ref. 14) through flow simulation of rough surfaces 
having Gaussian distribution of surface height. The 9x various 
roughness patterns is given by the following empirical relationship: 


= 1 - C e for Y < 1 


and 


(8) 


9x = 1 + C 


(r 


for 


> 1 


where C and r are constants and can be found in reference 14. The shear 
flow factor 95 is a function of h/a, the standard deviations, and the 
surface patterns of the surfaces a and b. Patir and Cheng (ref. 14) have 
shown that the shear flow factor 95 varies linearly with (ag/a)^ 
and (ob/a)^ and can be expressed in the following form: 



(9) 


The *s is given by 



exp 




and 


9 s = A 2 e 


-a^(h/a) 



for 


h 


71 


5 


5 


where l^i, A2, a^, 02, 03, and 04 are constants that can be found in refer- 
ence 14. 

The ratios (aa/a)*^ and (atj/a)^ are called the variance ratios of sur- 
faces a and b and are denoted by and Vpj^, respectively. The variance 

ratio of surface a can be written as 

Va - (^) = 1 - (^) (lO) 

since 



If both surfaces have the same surface pattern parameter, equation (9) 
reduces to 

■ <Va - ''rbl’s (?• i) 


By combining equations (10) and (11), q»s expressed as 

= (^Va - (?’ •') (12) 

When there is sliding, the effect of ic comes into the picture and 
becomes a parameter. If one of the surfaces is smooth, *5 is either a pos- 
itive or a negative quantity. The physical significance of $5 becomes 
clearer when the significance of the three values of variance ratio, namely, 
Vf-a = 1, 0.5, and 0, is understood. The variance ratio V,-a = 1 corresponds 
to the case shown in figure 2(a), where the rough surface is moving faster 



smooth one, V^g = 1; Ug > faster than rough one, 

$5 = +ve. Vpg = 0; Ug > U[j; = -ve. 



(c)Two surfaces have equal roughness. V^g =0.5; 
Ua ■ Ut)i 4>s " 0- 


Figure 2. - Surface roughness configurations in sliding. 
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than the smooth one. Conversely, V^a = 0 represents a case where the smooth 
surface is moving faster than the rough one (fig. 2(b)). And when = 0.5, 
the two surfaces have equal roughnesses (i.e., Og = 0 |j), *5 “ ^ (fig* 2 (c)).. 
Introducing equation (12) in to equation ( 6 ), we get 


3x 3x J 


3 h 3^ 

= 12“ T3T " 12“s»(2Va ‘ D a3T 


(13) 


It is shown in reference 1 that the expression for ahj/ax can be written 
in terms of nominal film thickness h as 


Hi. 1 

ax 7 


1 + erf 


fe) 


ah 

ax 


(14) 


Substitution of ahj/ax into equation (13) yields 


i_ , jL ia . i2u 1 

ax ^x 9x 2 


1 + erf 



a$ 

. 12u3,(2»„ - 1) 


(15) 


By using the substitutions X = x/R, H = h/R, Q = q/E‘, a = hQ/a, and 
Hq = hQ/R, equation (15) is nondimensionalized. The dimensionless form 
of equation (15) is written as 


aXlJ’x axy “ 
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FR I 


1 + erf 1 — 

,V2 H, 


IP /*^0\ 

lx (2V^g - 1) -air 


(16) 


To express in terms of no» the effect of temperature on viscosity 

has to be considered. 

Let us now put e = (aT/bq) into equation (5). The expression for 
ni then can be written as 


m = nn exp 


(1 


(17) 


where 


lb = 




Substituting equation (17) into equation (16) yields 



9 

3X 



3 3Q\ ^ 

aX/ " 2 
/ 


92 


1,3 aQ ^^1 
" aX aX 




aH 



( 18 ) 


where 


TIqU 

U = A 


and 


R = — 
s u 


Equation (18) can also be written in the following form: 

r2 


aX axj 


exp 


-■^b ll + 176 - 


T 


(1 - V) 


b m 3 aQ ae 

iXaX 


12u^ |- 


1 + erf 


aH 


\y/7 Hf 


aH / ^o\ 

ll * «sbr)(2Va - ^ 


b I 1 + T^e 


- 1 ' 


(19) 


To know H, we have to find the elastic deformation of the cylinders. 
The deformation v is calculated for a semi-infinite solid in a state of 
plane strain. The equation of deformation for two cylinders is (ref. 1) 


V(X) = - 


f ^2 

\ J ^ 

IT •'c 


P^)ln|X - S| dS 


(20) 


where 


V = 


IT 


8 



and 


S 


s 


= IT 


The procedure for calculating V is shown in detail in reference 1. After 
V is introduced in equation (3), the dimensionless film thickness H will 
be given by 


H = Hq + ^ + V (21) 


Before finding a solution of equation (19) for particular values of 
A, Hq, U, Rg, and Th, we must express e in terms of X. Using the modi- 
fied form of Blok's flash temperature model (ref. 7), we can express the tem- 
perature rise e (= aT/3q) by the equation 


9(X) 



Vx - s 


dS 


( 22 ) 


where 




and Pq is the dimensionless contact pressure, = Pr/E'- Therefore 0 can 
be obtained after evaluating the contact pressure. The contact pressure has 
been estimated by using a method given by Greenwood and Tripp (ref. 15) and 
by reference 1. 

Now knowing H, 9x» ^nd e, we can solve equation (14) by finite differ- 
ence methods for given values of G, U, R^, A, Hq, Tk, and T^.. The boundary 
conditions used are Reynolds boundary conditions. The details of the method 
can be found in reference 1. The Q distribution is transformed to P dis- 
tribution from 

P = - I In (1 - GQ) (23) 


Having obtained the pressure distribution, we can calculate the hydrodynamic 
load from 

^2 

P dx (24) 



or 


9 



( 25 ) 


= - I / ^ 1n(l - GQ) dX 
where 

Wu 

W - ^ 

" E'R 


The contact load is 



where 



and 


(26) 


(27) 


W 


c 


w 


c 


■ E'R 


The integrations of equations (25) and (27) are performed numerically by 
Simpson's rule. The total dimensionless EHL load Wj is 


«T 

" E'R 


(28) 


where 

Wt = w„ + w. 
T c h 


RESULTS AND DISCUSSION 

From the theory it can be seen that there are two kinds of dimensionless 
parameters - parameters relating to EHL and those due to the surface rough- 
ness effect including frictional heating. The EHL parameters are Wh, U, H, 
R^, and G. The other parameters are W^, a, a/6, y, T(-, and T^,. 

Thus both film thickness and temperature are functions of Wj, U, Rg, G, a, 


10 


Vra» ^b* the foregoing theory EHL load, film thickness 

and temperature rise were calculated for various speeds, slide-roll ratios, 
roughness parameters, variance ratios, surface patterns, and heating and ma- 
terials parameters and are presented in figures 3 to 12 and tables I to III. 
The parameters 6, U, Rs, a, a/p, and cover a wide range of practi- 

cal situations. In the following paragraphs the effect of some of these pa- 
rameters is discussed. 


Effect of Hydrodynamic Roughness Parameter 

The total EHL load decreased with increases in the surface roughness 
(fig. 3). (The smaller the values of A, the higher is the surface rough- 
ness). The minimum film thickness decreased slightly as the surface rough- 
ness was increased (fig. 4), and the maximum temperature rise increased with 
increases in the surface roughness (fig. 5). The variation of EHL load, min- 
imum film thickness, and maximum temperature with the roughness parameter a 
for various materials parameters G is also shown in figures 3 to 5 for cer- 
tain input conditions. The trend of EHL load and minimum film thickness with 
U, 6, and A was found to be similar to that observed when the temperature 
effect was not considered (ref. 1). In figure 6 the distribution of temper- 
ature rise is plotted for two roughness parameters (a = 1 and 2) and two ma- 
terials parameters (6 = 5000 and 7500). The maximum temperature rise always 
occurred in the inlet zone of the Hertzian contact. Similar observations have 
been made experimentally by Hamilton and Moore (ref. 9). Although in refer- 
ence 11 it was given that the maximum temperature rise would occur nearly in 
the exit zone, Jaeger (ref. 3) showed that his model could give the minimum 
temperature rise in the inlet zone. 


Effect of Variance Ratio 

The EHL load, minimum film thickness, and maximum temperature rise are 
plotted in figures 7 to 9 for three variance ratios (Vv-^ = 1, 0.5, and 0). 
The EHL load did not change significantly when V^^ changed from 0 to 0.5. 
But when the rough surface was moving faster than the smooth one (i.e., when 
Vpa = 1.0), the EHL load dropped (fig. 7). However, there was little or 
no variation in minimum film thickness with the change of variance ratio 
(fig. 8). The temperature rise was highest for = 1.0, lowest for 
Vp 3 = 0.5, and intermediate for = 0 (fig. 9). 


Effect of Surface Pattern Parameter y 

In figures 10 to 12 the variations in EHL load, minimum film thickness, 
and maximum temperature rise are shown for various surface pattern parame- 
ters. The surface pattern parameters y of 1/6, 1, and 6 represent trans- 
versely oriented, isotropic, and longitudinally oriented surfaces, respec- 
tively. For a particular set of input conditions, longitudinally and trans- 
versely oriented surfaces gave maximum and minimum EHL load, respectively 
(fig. 10). The smallest value of minimum film thickness was obtained for 
isotropic surfaces (fig. 11). The maximum temperature rise increased at a 
steeper rate with increasing a for the longitudinally oriented surface 
(fig. 12). The curves representing y of 1/6 and 1 are, for all practi- 
cal purposes, identical. 
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Materials 

Central 

parameter. 

film 

G 

thickness. 


Ho 

2500 

15. 0x10'* 

5000 

2as 

7500 

27.0 


xT 

o 


3 4 5 

Roughness parameter. A 


J 

6 



Figure X - Variation of EHL load with roughness parameter 
for various materials parameters, y = 1; o/p = Q. 0001; 


u = lO'll m/s; Rj - 0. 5; ■ 3000; T,, - 5; ■ 0. 5. 


6x10'^ 



Figure 5. - Variation of maximum temperature rise with 
roughness parameter for various materials parameters. 
r = l; a/p = 0.0001; U-lO'llm/s; Rs = 0.5; Tc = 3000; 
Tb'5i Vra-0.5. 


Materials 

Central film 

parameter. 

thickness. 

G 

Ho 

2500 

15.0x10"* 

5000 

2 a 8 

7500 

27.0 


30x10"* 



Figure 4. - Variation of minimum film thickness with rough- 
ness parameter for various materials parameters, y = 1; 
o/p = a 0001; U - lO'll m/S; R, ■ aS; T- ■ 3000; Th - 5; 

Vra-<^5. 
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10 ' 


rt 


Velocity. U; 
central film 
thickness, Hq 

10’’; 500x10-6 


10‘*0; 102.2xl0‘® 


Variance 

ratio, 

Vra 


1.0 



0 1.5 2.0 2.5 3.0 


Roughness parameter. A 

Figure 8. - Variation of minimum film thickness with roughness 
parameter for various variance ratios, y = 1; G = 5000; 
o/p = a0001; R5 = 0.5; 7^ = 3000; 7^ = 5. 


Maximum temperature rise. 



Roughness parameter, A 


Figure 9. - Variation of maximum temperature rise with roughness para 
meter for various variance ratios, y = 1; G = 5000; o/p = a 0001- 
Rs = a5; 7j. = 3000; 7u = 5. 





Figure la - Variation of EHL load with roughness parameter for various Figure a - Variation of minimum film thickness with roughness para- 

surface patterns. G = 5000; o/p = 0. 0001; Rg = 0. 5; = 3000; V^g “ OL 5. meter for various surface patterns. G = 5000; o/p = 01 0001- R,-CL5- 

T^= 3000; V„ = a5. 



Figure 12. - Variation of maximum temperature rise with rough- 
ness parameter for various surface patterns. G = 5000; 
a/p - aOOOl; Rj - a5; Tj. = 3000; \ = 5; Vra = aS. 
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Effect of Heating Parameter T^- 


In table I the variation of minimum film thickness and maximum tempera- 
ture rise for three heating parameters (Tc = 50, 3000, and 6000) is shown. 

- Although the minimum film thickness did not change appreciably with Tc, the 
temperatures rose with increases in the heating parameters. This is evident 
from the physical point of view. 


Effect of Surface Contact Roughness Parameter a/3 

The higher the value of a/3, the rougher is the surface. The rise in 
maximum temperature increased with increases in the value of a/ 3. The 
minimum film thickness also fell slightly as a/3 was increased (table II). 


Effect of Slide-Roll Ratio Rg 

The slide-roll ratio Rs = 0 corresponds to a pure-rolling case. The 
effect of Rg on and for various is given in table III. 

When both surfaces had the same roughness structure (i.e., = 0.5), both 

minimum film thickness and maximum temperature rise did not change with the 
slide-roll ratio. The minimum film thickness decreased as the slide-roll 
ratio was increased for variance ratios V^a of 0 and 1.0. On the other 
hand the maximum temperature rise dropped when = 0 and increased when 

^ra “ 1*0. 


EXAMPLE 

From the temperature distribution 0 it is sometimes difficult to es- 
timate the quantitative value of temperature rise in the contact zone. The 
following example illustrates a calculation procedure for finding aT. 

Two heavily loaded cylinders having surface a rough and surface b 
smooth are lubricated under the following conditions: 

Cylinder radii, Rg = 0.033 m, Rj^ = 0.05 m 
Velocity of surface of cylinder a, Ug = 9 m/s 
Velocity of surface of cylinder b, U 5 = 4 m/s 

Mean coefficient of friction for a contact between asperities, u = 0.5 
Equivalent modulus of elasticity of cylinder materials, 

E' = 2.2xl0ll N/m2 

Density of material, p = 8x10^ kg/m^ 

Thermal conductivity of material, k = 10 J/m K s 
Specific heat of material, c = 100 J/kg K 
Central film thickness, ho = 0.416x10“® m 

Pressure-viscosity coefficient of lubricant, a = 2.27x10“° m^/N 
Viscosity of lubricant, no = 0.0068 Pa s 

From these data and 30 = 2050 K, the following parameters are calculated: 

Dimensionless speed parameter. 

U s 10“11 
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TABLE I. - VARIATION OF MINIMUM FILM THICKNESS AND MAXIMUM TEMPERATURE 
RISE FOR VARIOUS HEATING PARAMETERS 
[y = 1; G = 5000; o/3 = 0.0001; Rg = 0.5; T^ = 5; V^^ = O.5.] 


Surface 

velocity, 

u. 

m/s 

Cantral 

film 

thickness, 

Ho> 

m 

Hydro- 

dynamic 

roughness 

parameter, 

A 


Heating parameter 

Tc 



50 

3000 

6000 

Values of minimum film thickness, m, and maximum 
temperature rise, K 

^min 

®max 

^min 

®max 

^min 

®max 

10x10-12 

20.80x10-6 

1 

17.96x10-6 

0.438x10-3 

18.39x10-6 

2.597x10-3 

17.87x10-6 

5.440x10-3 



2 

17.76 

.003 

17.79 

.220 

17.76 

.398 



3 

20.23 

0 

20.23 

.001 

20.23 

.003 


TABLE II. - VARIATION OF MINIMUM FILM THICKNESS AND MAXIMUM TEMPERATURE RISE FOR VARIOUS RATIOS 
OF PRESSURE-VISCOSITY COEFFICIENT TO MEAN RADIUS OF CURVATURE OF ASPERITIES 
[y = 1; G = 5000; Rg = 0.5; T^. = 3000; T^ = 5.] 


Surface 

velocity, 

u, 

m/s 

Central 

film 

thickness, 

Ho. 

m 

Hydro- 

dynamic 

roughness 

parameter, 

A 

1—1 

o 

o 

o 

o 

II 

CO. 

D 

o/B = 

0.01 

Values of minimum film thickness, m, 
maximum temperature rise, K 

and 

^min 

®max 

^min 

®max 

10x10-12 

20.80x10-6 

1 

18.39x10-6 

2.289x10-3 

18.35x10-6 

3.412x10-3 



2 

17.79 

.220 

17.72 

.224 



3 

20.23 

.001 

20.39 

.014 












TABLE III. - VARIATION OF MINIMUM FILM THICKNESS AND MAXIMUM TEMPERATURE 
FOR VARIOUS SLIDE-ROLL RATIOS AND VARIANCE RATIOS 
[y = 1; G = 5000; a/e = 0.0001; U = 10-H m/s; Hq = 20.8xl0"6 m; Tc = 3000, Tb = 5.] 


Variance 

ratio, 

Vra 

Hydro- 

dynamic 

roughness 

parameter, 

A 


Slide-roll 

ratio, Rs 


0 

0.5 


Values of minimum*film thickness, 
and maximum temperature rise, K 

n. 

^min 

®max 

^min 

®max 

0 

1 

18.40x10-6 

2.397x10-3 

17.48x10-6 

2.371x10-3 


2 

17.79 

.223 

17.55 

.223 


3 

20.23 

.001 

20.25 

.001 

.5 

1 

18.40 

2.289 

18.40 

2.289 


2 

17.79 

.218 

17.79 

.218 


3 

20.23 

.001 

20.23 

.001 

1.0 

1 

18.40 

2.289 

18.13 

3.231 


2 

17.79 

.218 

18.59 

.169 


3 

20.23 

.001 

20.11 

.002 


Dimensionless central film thickness. 

Hq = 20.8x10-6 

Dimensionless materials parameter. 

G £ 5000 

Heating parameter. 

Tc *5? 3000 

Assuming y = 1» A = 1, Rg = 0.5, Tb = 5, and o/b = 0.0001, vt/e get from 
figure 9, = 3.2x10“^. Now the maximum temperature rise AT^ax = 

®max^0 6.5 d K. 


CONCLUSIONS 

From the data presented herein the following conclusions were drawn: 

1. The elastohydrodynamic load for a constant central film thickness 
decreases with increasing hydrodynamic surface roughness. 

2. Minimum film thickness decreases with increasing roughness. 

3. The maximum rise in temperature increases with increasing surface 
roughness. 

4. The location of maximum temperature is always in the inlet zone. 
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5. The minimum film thickness decreases as the cylinders start sliding from 

pure rolling. ^ ^ u n • ^ n 

6. Film thickness and temperature rise are independent of the slide-roll 

ratio for surfaces with equal roughness structures. 


7. When the rough surface is moving faster than the smooth one (variance 
10 Vy,a = 1 ), the temperature increases for higher slide-roll ratios. 

8. There is no significant change in minimum film thickness when 


rati 


varies from 0 to 1. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 2, 1981 
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